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Abstract 

A rigorous quantum relativistic approach has been used to calculate the re- 
lationship between the decay laws of an unstable particle seen from two inertial 
frames moving with respect to each other. In agreement with experiment, it 
is found that the usual Einstein's time dilation formula is rather accurate in 
this case. However, small corrections to this formula were also obtained. Al- 
though the observation of these corrections is beyond the resolution of modern 
experiments, their presence indicates that special relativistic time dilation is not 
rigorously applicable to particle decays. 

1 Introduction 

Unstable particles are a perfect testing ground for theories that try to unify the princi- 
ple of relativity with quantum mechanics for (at least) two reasons. First, an unstable 
particle provides a very simple example of a non-trivial interacting quantum system. 
Due to the specific initial condition characteristic to this problem (there is just one 
particle in the initial state), a rigorous description of the decay is possible in a small 
Hilbert space that contains only states of the particle and its decay products, so the 
solution can be obtained in a closed form. Second, unstable particles played an impor- 
tant role in confirming predictions of Einstein's special relativity. It was demonstrated 
experimentally that the decay of moving particles slows down in a good agreement 
with Einstein's time dilation formula PjJ |2] . 

For a long time it was believed that relativistic quantum mechanics and quantum 
field theory must reproduce exactly the Einstein's time dilation formula for unstable 
particles [Sj. However, a detailed quantum relativistic calculation in [I] has shown 
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that there are small corrections to this formula in the case of particles with definite 
momentum. This result was later confirmed in other studies jSJ Ej. Although sur- 
prising, these findings did not challenge directly the applicability of special relativity 
to unstable systems. This is because relativistic transformations connect results of 
measurements of two inertial observers moving with constant velocities with respect to 
each other. Then a fair comparison with the time dilation formula requires considera- 
tion of unstable states having definite values of velocity for both observers. Such states 
do not have a definite momentum for, at least, one observer. Their decay laws were 
considered in recent article [Zj which alleges that the decay accelerates when particle's 
velocity increases, instead of being slowed down as experiment shows. In the present 
article we resolve this controversy by performing a detailed calculation of the decay 
laws of particles with narrow distributions of velocities observed from moving reference 
frames. We found that the result of [7] is based on an incorrect identification of the 
subspace of states of the unstable particle in the full Hilbert space of the system. In 
sections - 13 we present a rigorous quantum relativistic framework required for the 
description of decays. The exact formula for the time dependence of the non-decay 
probability in a moving inertial frame of reference is derived in section |HJ Particular 
cases of this formula relevant to unstable particles with sharply defined momenta or 
velocities are considered in sections El and [TTJ respectively. Section H21 is devoted to 
numerical calculations of the differences between the accurate quantum mechanical 
result and the standard time dilation formula ©. Although these differences are much 
smaller than the resolution of modern experiments, their presence is in sharp contra- 
diction with special relativity which does not tolerate even small deviations from the 
Einstein's time dilation. 

2 Formulation of the problem 

The decay of unstable particles is described mathematically by the non-decay proba- 
bility which has the following definition. Suppose that we have a piece of radioactive 
material with N unstable nuclei prepared simultaneously at time t = and denote 
N u (t) the number of nuclei that remain undecayed at time t > 0. Then the non- 
decay probability 1 u(0,t) (also called the non-decay law in this paper) is defined as 
the fraction of nuclei that survived the decay in the limit of large N 

u(0,t) = lim N u (t)/N (1) 

iV— »oo 

So, at each time point the piece of radioactive material can be characterized by its 
composition (N u (t)/N; (1 — N u (t)/N), where N u (t)/N is the share of undecayed nuclei 

1 The first argument in u(0,t) indicates the rapidity of the observer that measures the non-decay 
probability, as described below. 
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in the sample and 1 — N u (t)/N is the share of nuclei that were transformed to the 
decay products. In this paper, in the spirit of quantum mechanics, we will treat iV 
unstable particles as an ensemble of identically prepared systems and consider u(0,t) 
as a property of a single particle (nucleus), i.e., the probability of finding this particle 
in the undecayed state. 

In this paper we are concerned with comparison of decay observations made by two 
observers O and O' that move with respect to each other. Without loss of generality 
we will assume that observer O' moves with respect to O with velocity (0,0, v) along 
the z-axis and that at time t — measured by both observer's clocks the origins of 
their coordinate systems coincide and all three pairs of coordinate axes are parallel 
(x || x', etc.). For discussion purposes we will say that O is at rest while O' is moving. 
To simplify formulas we will use the rapidity parameter 6 = tanh _1 (t>/c) instead of 
velocity v. For example, in this notation the famous relativistic factor takes the form 
7 = (1 — t> 2 /c 2 ) _1 / 2 = cosh#. If from the point of view of O the decay is described by 
the function u(0,t), then the non-decay law of the same particle seen by the moving 
observer O' will be denoted by u>(9, t). 2 Calculation of this function and, in particular, 
the relationship between u(9,t) and u(0,t) is the major goal of this work. 

3 Postulates of relativity 

The purpose of any relativistic theory is to describe the relationships between measure- 
ments made by different inertial observers or from different inertial reference frames, 
i.e., reference frames moving in space with constant velocities along straight lines and 
without rotation. The special principle of relativity tells us that these reference frames 
are exactly equivalent 

• Postulate I. Experiments identically arranged and performed in two different 
inertial reference frames O and Q always yield the same results. 

For any two observers O and Q there is an inertial transformation that connects Q to 
O, i.e., the set of rules that allows one to change from the reference frame O to the 
reference frame Q. Each inertial transformation is a combination of space and time 
translations, rotations and boosts. A composition of two inertial transformations is 
again a valid inertial transformation. This composition obeys the associativity law. 
The inverse of an inertial transformation is also an inertial transformation. Therefore, 
they form a 10-parameter Lie group. The structure of this group is specified in the 
second postulate 

• Postulate II. Inertial transformations (space and time translations, rotations, 
and boosts) form the Poincare group. 

2 Note that here parameter t denotes time measured by the clock belonging to the observer O'. 
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Obviously, two different observers obtain different results by measuring observables 
of the same physical system. Most problems in physics can be understood as trans- 
lations of descriptions of the physical system between different reference frames. For 
example, if we have a full description of the system from the point of view of observer 
O, then the time evolution is obtained by answering the question "what is the descrip- 
tion of the same system from the point of view of observer O" that is shifted in time 
with respect to O?" It is important to realize that inertial transformations between 
observers can be divided into two groups: kinematical and dynamical. Kinematical 
transformations are those whose action on observables does not depend on the inter- 
action. For example, if two stationary observers look at the unstable particle from 
different points in space, they would assign the same non-decay probability to the par- 
ticle. The same is true for stationary observers having different orientations in space. 
From this follows 3 

• Postulate III. Space translations and rotations are kinematical. 

On the other hand, time translations produce non-trivial changes in the system. The 
composition of the unstable system looks different for observers O and O" shifted in 
time with respect to each other. The exact action of time translations on observables 
of the physical system should be obtained as a result of solution of dynamical equations 
which depend on the interaction acting in the system. Hence the following postulate 
is true for all types of isolated interacting systems. 

• Postulate IV. Time translations are dynamical, i.e., interaction-dependent. 

The above postulates I - IV have overwhelming experimental support. They are 
assumed to be valid throughout this paper. However they are not sufficient for a full 
description of the unstable particle in a moving reference frame. Such a description 
requires also knowledge of the nature of boosts. Here we have a choice between two 
paths forward. One path is to postulate certain properties of boosts. This path was 
taken by Einstein. It leads to special relativity and to the time dilation formula (jHJ) as 
explained in the next section. In this paper we will argue in favor of choosing another 
path: keep postulates I - IV, add to them well-established postulates of quantum 
mechanics, and see what are the implications for the transformations of observables 
(in particular, the non-decay probability) with respect to boosts. This approach is 
employed starting from section El throughout the paper. 

3 Our adoption of this postulate means that we are working in Dirac's instant form of relativistic 
dynamics Other forms of dynamics, e.g., the point form and the front form, were also introduced 
by Dirac and are frequently used for description of relativistic interactions. As discussed in section 
1131 these forms are not appropriate for the description of non-decay laws of unstable particles. 
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4 Particle decay in special relativity 



In addition to postulates I - IV, Einstein's special relativity makes two more assump- 
tions regarding the nature of boosts 



• Assumption V. Boosts are kinematical. 



By postulating the kinematical character of boosts, special relativity insists that the 
internal composition of a compound system does not depend on the velocity of the 
observer. For example, if O and O' are two observers moving with respect to each other, 
then, according to special relativity, both observer will measure the same composition 
of the unstable system at time t = 0. These statements are often considered as self- 
evident in discussions of special relativity. For example, R. Polishchuk writes in [5| 
"Any event that is "seen" in one inertial system is "seen" in all others. For example if 
observer in one system "sees" an explosion on a rocket then so do all other observers. " 
In addition to the universality and interaction-independence of boost transformations, 
special relativity also postulates the exact transformation laws of physical observables 
with respect to boosts. They are referred to as Lorentz transformations. The most 
fundamental are Lorentz transformations for space-time coordinates of events. 4 

• Assumption VI. If from the point of view of observer O an event is localized 
in a space point (x, y, z) at time t, then from the point of view of observer O' the 
same event has space-time coordinates (t f , x', y', z') given by Lorentz formulas 



t' = t cosh 9-- sinh 9, (2) 

c 



x' = x, (3) 

y' = y, (4) 

z = z cosh 6* — ctsinh^, (5) 



Let us now demonstrate that in special relativity the Postulates I - IV and As- 
sumptions V - VI are sufficient to unambiguously describe particle decay in different 
reference frames without invoking any information about the interaction governing the 
decay. Suppose that from the point of observer O the unstable system is prepared in 

4 By event we understand a measurable physical process occurring at a certain point in space at 
one time instant. An intersection of trajectories of two point-like classical particles is an example of 
such an event. Note that in quantum mechanics the definition of event is problematic as particles do 
not have well-defined trajectories. Moreover, it is known that a particle sharply localized from the 
point of view of the observer O loses its localization from the point of view of the moving observer O' 
So, Assumption VI can be applied only in the classical limit. 
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the state with composition (1.0; 0.0) at rest in the origin x = y = z = 0a.t time t = 0. 5 
Then observer O may associate the space-time point 



(t, x, y, z) pre p = (0, 0, 0, 0) (6) 

with the event of preparation. In accordance with the dynamical character of time 
translations, the non-decay probability u(0,t) decreases with time. From experiment 
and quantum mechanical calculations (see, e.g., section [T2"|) it is known that the non- 
decay law has an (almost) exponential shape 6 



u(Q,t) « exp( ) (7) 

TO 

where tq is the lifetime of the unstable particle. At time t — To the non-decay prob- 
ability is exactly e" 1 , so that the composition is (0.368; 0.632). This "one lifetime" 
event has coordinates 



(t,x,y,z) Hfe = (t , 0,0,0) (8) 

according to the observer O. 

Let us now take the point of view of the moving observer O'. Due to the Assumption 
V, both O and O' agree that at the "preparation" and "lifetime" events the non-decay 
probabilities are 1 and e _1 , respectively. However, observer O' may not agree with 
O about the space-time coordinates of these events. Substituting (0) and (JBJ) in (J2J) 
- © we see that from the point of view of O', the "preparation" event has coordi- 
nates (0,0,0,0), and the "lifetime" event has coordinates (r cosh^, 0, 0, — cro sinh 9). 
Therefore, the time elapsed between these two events is cosh 9 times longer than in 
the reference frame O. This also means that the decay law of the particle is exactly 
cosh 9 slower from the point of view of the moving observer O' . This is reflected in the 
famous "time dilation" formula 



This formula was confirmed in numerous experiments []], most accurately for muons 
accelerated to relativistic speeds in a cyclotron These experimental findings were 

5 As indicated in the previous footnote, here we use the classical limit. In quantum mechanics the 
notions "at rest" and "in the origin" cannot be used simultaneously. 
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The exact form of the non-decay law is not important for our derivation of eq. here. 
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certainly a triumph of Einstein's theory. However, as we see from the above discussion, 
in special relativity eq. Q can be derived only under two assumptions V and VI, 
which lack proper justification. Therefore, a question remains whether eq. (|9|) is a 
fundamental exact result or simply an approximation that can be disproved by more 
accurate measurements? 

Our goal in the rest of this paper is to demonstrate that this classical result of 
special relativity is not exact. In section H21 we will calculate corrections to the formula 



5 Quantum mechanics of particle decays 

Let us now turn to the description of an isolated unstable system from the point of 
view of relativistic quantum theory. In our approach we will keep postulates I - IV. 
However, we are not going to use Assumptions V and VI. 

We will consider a model theory with particles a, b, and c, so that particle a is 
massive spinless and unstable, while its decay products b and c are stable and their 
masses satisfy the inequality 



which makes the decay a — *> b+c energetically possible. In order to simplify calculations 
and avoid being distracted by issues that are not relevant to the problem at hand we 
neglect the spin of the particle a and assume that there is only one decay mode of this 
particle, i.e., into two decay products b and c. For our discussion, the nature of the 
particle a is not that important. For example, this could be a muon or a radioactive 
nucleus or an atom in an excited state. 

When measuring the non-decay law, experimentalists simply count the number of 
particles and determine their types. For example, observations of a muon may result in 
only two outcomes. One can find either a non-decayed muon or its decay products (an 
electron, a neutrino, and an antineutrino). Thus the number of particles is a legitimate 
observable, and for our model system a — > b + c we can introduce Hermitian operators 
for the number of particles N a , Nb, and N c . Since these observables can be measured 
simultaneously, we can assume that these three operators commute with each other. 
Two combinations of their common eigenvalues are allowed in our system: 



m a > m.b + m c 



(10) 



n a = l,n b = Q,n c = 
n a = 0,n 6 = l,n c = 1 




(12) 
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Hence the Hilbert space of the unstable system should be represented as a direct sum 
of two orthogonal subspaces 7 



Ti = H a ® Hbc 



(13) 



where 7i a is the subspace of states of the unstable particle a which corresponds to the 
set of eigenvalues (JUJ), and 7ib c = Hb <S) TC C is the orthogonal subspace of the decay 
products which corresponds to the set of eigenvalues (fT2j) . 

We can now introduce a Hermitian operator T (also known as "yes-no experiment" ) 
that corresponds to the observable "particle a exists". The operator T can be fully 
defined by its eigensubspaces and eigenvalues. When a measurement performed on the 
unstable system finds it in a state corresponding to the particle a (i.e., the state vector 
is within Ti. a ), the value of T is 1. When the decay products b + c are observed (the 
state vector lies in TCb c ), the value of T is 0. Apparently, T is a projection operator on 
the subspace Ti. a . For each normalized state vector GW the probability of finding 
the unstable particle a is given by the expectation value of the observable T [3] 



In relativistic quantum mechanics, the dynamics of the system is described by 
a unitary representation U g of the Poincare group in the Hilbert space H. [TT] . It is 
convenient to express representatives U g as exponential functions of generators P, J, K, 
and H. Then a general unitary operator from the set U g can be always written as a 
product 



where <p is the rotation vector, v = tanh^ is the velocity of the boost, a is the vector 
of space translation, and t is the amount of time translation. Generators are Hermitian 
operators and correspond to certain total observables of the system. Space translations 
are generated by the vector of the total linear momentum P. The generator of rotations 
J is the operator of the total angular momentum. The generator of time translations 
H is the Hamiltonian (total energy). The generator of boosts K is called the boost 
operator. These generators must satisfy the commutation relations of the Poincare 
group Lie algebra. In this paper we will need, in particular, the following commutators 

7 In principle, a full description of systems involving these three types of particles must be formu- 
lated in the Fock space where integer eigenvalues n a ,rib, and n c are allowed to take any values from 
zero to infinity. However, for most unstable particles the interaction between the decay products in 
the final state can be ignored, and considering the subspace l|13(l of the full Fock space is a reasonable 
approximation. 




(14) 
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[K i ,Pj\ = -i^HS ij (15) 

[K i ,H] = -ihP i (16) 
where i,j = x, y, z. They imply the following useful relationships 8 



e ^ Pe ^Kce = p + £[(p. ^)( C osh0-l) + ~Fsinh0] (17) 

9 9 c 

e i Kci He-i^ = H coshfl + c(P - B -) sinhfl (18) 

The representation (fTlj) allows us to relate results of measurements in different 
reference frames. Let us first take the point of view of the observer O and consider a 
vector = | \I/ (0, 0) ) G TC a 9 that describes a state in which the unstable particle a is 
found with 100% certainty. 



^(0,0) = (tt|T|tt) 
= 1 

Then the time evolution of the state vector in the reference frame O is given by 



|tt(0,t)) = e-* H *|*> 
and the non-decay law is given by 



u m (0,t) = (V\e£ m Te-£ m \y) 

From this equation it is clear that the Hamiltonian H describing the unstable system 
should not commute with the projection T 



[H,T}^0 (19) 

Otherwise, the subspace 7i a of states of the particle a would be invariant with respect 
to time translations and the particle a would be stable (o>|^)(0, t) = 1 for all t). 

8 For derivation see section 2.2 in IT21 . 

9 In the notation 1^(0, 0)) the first argument is the rapidity parameter of the reference frame from 
which the state |\&) is observed and the second argument is the time of observation. This is consistent 
with the convention adopted for uj(0,t) in section^ 



9 



The moving observer O' describes the initial state (at t = 0) by the vector 

|tf(0, 0)) = e^ e \^>) 
The time dependence of this state is 

|tf(0,t)) = e-^ m e^ e \^>) (20) 
Then the non-decay law from the point of view of O' is 

u m (6,t) = (*(6,t)\T\*(8,t)) (21) 

= \\Tme,t))\\ 2 (22) 

where the last equation follows from the property T 2 = T of the operator T. 

6 Non-interacting representation of the Poincare 
group 

Before calculating the non-decay law in the moving reference frame (|21|) . let us first 
consider a simpler case when the interaction responsible for the decay is "turned off" . 
Then the representation of the Poincare group U® acting in TC is non-interacting. This 
representation is constructed in accordance with the structure of the Hilbert space (|T3|l 
as 

U° g = U; © U b g ® U c g (23) 

where U^ a ' b ' c ^ are unitary irreducible representations of the Poincare group correspond- 
ing to particles a, b, and c, respectively. The generators of the representation (J23j) are 
denoted by P , Jo, H Q , and K . The operator of non- interacting mass 

Mo = +±jH$-PZc2 

commutes with P , Jo, H , and K . According to (fTUj). the operator M has a contin- 
uous spectrum in the interval [mj + m c , oo) and a discrete point m a embedded in this 
interval. 
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From definition ()23|) it is clear that the subspaces 7i a and 7ih c are invariant with 
respect to C/°. Moreover, the particle number operator N a and the projection operator 
T commute with the non-interacting generators 

[T, P ] = [T, J ] = [T, K ] = [T, H ] = (24) 

This implies that the particle a is stable with respect to time translations and boosts, 
as expected 

u m (6,t) = (^\e-f^ e e^ Hot Te-^ Hot ef (K ^ e \^) 
= ($\T\9) 
= 1 

The primary reason for considering the representation Ug is that it allows us to 
build a convenient basis in the subspace 7i a . Let us denote |0) a vector in 7i a that 
corresponds to the particle a with zero momentum. 



Po|0) = 
H \0) = m a c 2 \0) 



Then we find that the vector 



|p) = e -n*° a \0) (25) 

describes the particle a with definite momentum p = |m a csinh^. Indeed, using (JT7J) 
we obtain 



Po|p) = P e-^ K ° e |0) 

= e -f K °^f K °0p oe -f K °^|o) 

= e-^ Ko0 \-p o + ^[(P ■ ^)(cosh^ — I) + -F sinh^])|0) 
c 



= -m a c sinh 0e~^ K ° e \ 0) 


Q 

= -m a c siiah 0\p) 


Since particle a is spinless by our assumption, the eigenvectors |p) of the total momen- 
tum operator P form a full basis in the subspace TC a , so that 
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(p|p') = *(P-P') (26) 
T = y dp|p)<p| (27) 

Then any state \^f) G 7i a of the particle a can be represented by a linear combination 
of these basis vectors 



= y dp|p><pi*> 

= /dp^(p)|p> (28) 

where V'(p) = (p|^) i s the wave function in the momentum representation. In order 
to ensure the normalization (^l^) = 1, the function ip(p) must satisfy 

/dp|^(p)| 2 = l (29) 

Apparently, vectors |p) themselves are not normalized. If we want to study the non- 
decay law of a state with definite momentum p , we should use a state vector (which 
we denote by |p ) to distinguish it from |p ) ) that has a normalized momentum-space 
wave function sharply localized near p . In order to satisfy eq. (PSj) such a wave 
function may be formally represented as a square root of the Dirac's delta function 

v>(p) = \A(p - po) 

The action of boosts on the basis vectors |p) and |p) is known from Wigner's theory 
of irreducible unitary representations of the Poincare group (see, e.g., [TT^ ITT]) 

e-fW|p> = /If I^P> ( 3 °) 

where 

Lp = (p x ,Py,Pz cosh^ H -sinh^) 

f2 p = \J m^c 4 + c 2 p 2 
Using eq. (JHUJ) and the property 
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dp = d(Lp) ^ 



we can find boost transformations for an arbitrary state vector of the form (J2%|) 

e -f<*.)««|*) = /^(pje-f^lp) 



/«(P)^|LP) 



These transformations can be viewed as transformations of the corresponding momentum- 
space wave function, e.g., 



e-f^V(p) - /^^(^p) ( 32 ) 

so that the boost operator can be represented as a differential operator in the momen- 
tum space 

W(P) = -fm^M^) 
c 0^0 cw 

= %^./& fc ,p^coshe-^sinhe) 
c e^o eztV y ii p c 

- -^+!t Wp) (33) 

For further calculations we will need to define the Newton- Wigner position operator 
PH in H 



which has the property 

[(Ri)o,(PM = ihSij (34) 
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that can be verified by direct substitution. According to eq. (J3~3~j) 



dp x 

Therefore operator e^^ R °^ zb acts as a translation operator in the momentum space. In 
particular, we can write 



e f (Bo Wsinh<y 5(p) = ^ (p _ po) 

where po = (0, 0, macsinh^). On the other hand, applying the boost transformation 
to the momentum eigenfunction we obtain 10 



'(J(p-Po) 

This suggests that momentum eigenvector has another useful representation 

|p) = e ^ R °- p |0) (35) 

7 Interacting representation of the Poincare group 

Let us now "turn on" the interaction responsible for the decay and discuss the inter- 
acting representation U g of the Poincare group in TC with generators P, J, K, and 
H. According to our postulates III and IV, the generators of space translations and 
rotations are interaction-free, 



10 where | J\ = — ^ lp is the Jacobian of the transformation p — > L 1 p. 
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P = Po 
J = Jo 



while the generator of time translations (the Hamiltonian H) contains an interaction- 
dependent term V. 



H = H Q + V 

From the commutator (|15p it then follows jHj that the generators of boosts must be 
interaction-dependent as well 



K = K + W 



where W/0. This means that we are working in the so-called Dirac's instant form of 
dynamics. We will further assume that the interacting representation U g belongs to the 
Bakamjian-Thomas form of dynamics |14j . which is characterized by the property that 

the interacting operator of mass M = c~ 2 \J H 2 — PqC 2 commutes with the Newton- 
Wigner position operator 11 



[Ro, M] = (36) 

Our next goal is to define the basis of common eigenvectors of commuting operators 
Po and M in Ti. 12 These eigenvectors must satisfy conditions 



P |p,m) = p|p,m) (37) 
M\p,m) = m\p,m) (38) 

They are also eigenvectors of the interacting Hamiltonian H — \l M 2 c 4 + Pf)C 2 

11 The possibilities for the interaction to be not in the Bakamjian-Thomas instant form are discussed 
in section 

12 In addition to these two operators, whose eigenvalues are used for labeling eigenvectors |p, m), 
there are other independent operators in the mutually commuting set containing Po and M . These 
are, for example, the operators of the square of the total angular momentum Jq and the projection 
of the total angular momentum on the z-axis (Jq) z - Therefore a unique characterization of any basis 
vector requires specification of all corresponding quantum numbers as |p, to, j 2 ,j z , . . .). However these 
quantum numbers are not relevant for our discussion and we omit them. 
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H\p,m) = u p \p,m) 

where u p = y/m 2 ^ + c 2 p 2 . 13 In the zero-momentum eigensubspace of the momentum 
operator P we can introduce a basis |0, m) of eigenvectors of the interacting mass M 

P |O,ra) = 
M\0,m) = m\0,m) 

Then the basis |p, m) in the entire Hilbert space 7i can be built by formula (cf. eq. 

|p, m) = e » \0,m) 
where p = mcQQ~ x sinh 9. These eigenvectors are normalized to delta functions 

(q, m\p,m') = <5(q — p)8(m — m) (39) 

The actions of inertial transformations on these states are well-known [TT|. In partic- 
ular, for boosts along the z-axis (cf. eq. (j30J) ) and time translations we obtain 



e-f K ' e \p,m) = J^\Ap,m) (40) 

V u p 

e^ Ht \p } m) = e^ pt \p,m) (41) 



where 



Ap = (Px,P v ,Pz cosh 6* H — -sinh^) 

Next we notice that due to eqs. ()34|) and ()36|) vectors e^ Ro ' p |0,m) also satisfy eigen- 
vector equations (|3*Tj) - (jSHJ), so they must be proportional to the basis vectors |p, to) 



|p, m) = 7(p, m)e n °' p |0,m) 

13 Note the difference between ui p that depends on the eigenvalue m of the interacting mass operator 
and Q p in eq. that depends on the fixed value of mass m a of the particle a. 
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Figure 1: Mass distribution of a typical unstable particle. 



where j(jp,m) is an unimodular factor. Unlike in (|H5|) . we cannot conclude that 
7(p,m) = 1, because, generally, the action of e^ R °' p on eigenvectors |q, m) involves 
multiplication by a unimodular scalar in addition to the shift of momentum. However, 
if the interaction is not pathological we can assume that the factor 7(p,m) is smooth, 
i.e., without rapid oscillations. This property will be used in derivation of eq. (JoTJj) . 

Obviously, vector |0) can be expressed as a linear combination of zero-momentum 
basis vectors |0,m), so we can write 14 



oo 

|0) = J dmn(m)\Q,m) (42) 

m b +m c 

where function |/i(m)| 2 describes the mass distribution of the unstable particle. For 
most realistic unstable systems the mass distribution \fi(m)\ 2 has the Breit-Wigner 
form (see Fig. I) 15 



2 aT/2Ti 

K m ~ r2 /A i 7 ^,tfm>m b + m c 43 

L z /A+ [m — vn a y 



14 We will assume that interaction responsible for the decay does not change the spectrum of mass. 
In particular, we will neglect the possibility of existence of bound states of particles b and c, i.e., 
discrete eigenvalues of M below mi, + m c . Then the spectrum of M (similar to the spectrum of Mq) 
is continuous in the interval [rrib + m Cl oo), and integration in l|42|) should be performed from mj, + m c 
to infinity. 

15 Strictly speaking, the center of the resonance 143(1 could be different from the mass m a that 
particle a has in the absence of decay interaction. However, we will disregard this possibility here. 
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\/i(m)\ 2 = 0, if m < m b + m c 



(44) 



where parameter a is a normalization factor required to ensure that the mass distri- 
bution (|4~H|) - (^1} is normalized to unity. 

We now use eqs. (|H5J) and ()42)1 to expand the vector |p) in the basis |p, m) 



|p) = e» RoP J dmfi(m)\0,m) 

m b +in c 

oo 

= J dmjj(m)'j(p,m)\p,m) (45) 

m b +m c 

Then, from (}39|) we obtain a useful formula 

oo 

(q|p, m) = J dm' n* (m')'j* (q, m')(q, m'\p, m) 

m b +m c 

= 7*(p, m)/i*(m)5(q — p) (46) 

8 General formula for the non- decay law 

Suppose that vector in (|2H|) describes a state of the unstable particle a from the 
point of view of the observer O. The time dependence of this state seen from the 
moving reference frame O' is obtained by applying eqs (J20J), (J4*U|) . (}4~T]) . and (|43j) 

oo 

= y dpijj(p) J dmji(m)^(p,m)e~^ Ht e^ Kze \p,m) 

m b +m c 

oo , 



= J dpip(p) J dm//(m)7(p, m)e lpt \j — 1 Pi m ) 

m b +m c ' 

The inner product of this vector with |q) is found by using (j4T?j) 
(q|*(0,t)> 

oo 

= J dpip(p) J dmfjL(m)j(p, m)e _ s a ' A_1 p t (q|A -1 p, m) ' " ' p 



m b +m c 
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= fdpip{p) J dm| / u(m)| 2 7(p,m) 7 *(A" 1 p,m)e-^- 1 p t 5(q- A- 1 ?)./^^ 

m b +m c ' p 

Introducing new integration variables r = A x p and taking into account this 
equation can be rewritten as 



<q|*(M)> 

./— ^(Ar)7(Ar)7*(r)|/i(m)|V^ r *(5(q - r) 

U r V ^Ar 



J dm 


[dr U - 


m b +m c 




oo 






/^Aq 


m b +m c ' 





The non-decay probability in the reference frame O' is then found by substituting (|27|) 
in eq. (JUJ) 

u>|*>(0,t) = |dq<*(0,t)|q)(q|*(0,t)) 
'dq|<q|¥(0,t)> | 2 



y dq| y rfm /^^(Aq)7(Aq,m)7*(q,m)|/i(m)|V^ q '| 2 (47) 



which is an exact formula valid for all values of 9 and t. 

9 Decay law in the reference frame at rest 

In the reference frame at rest (8 = 0), formula ()47|) simplifies 



oo 

1 2 1 



u\^)(0,t) = J dq\ip(q)\ \ J dm\n(m)\ 2 e * L 

m b +m c 

Let us consider a particular case of this expression when the unstable particle has a 
well-defined momentum, i.e., described by the state vector |p). If the particle a is in 



the state |0) with zero momentum then ?/>(q) = y<5(q) and from the above equation 
we obtain 



19 



W| O) (0,t) = | J dm|/i(m)| 2 e-> c2 f (48) 

m b +m c 

which is the standard formula for the non-decay law of a particle at rest. 16 In particular, 
if we substitute fl4*Hj) - (jH)) for the mass distribution |/x(m)| 2 we obtain approximate 
exponential decay (J7j) where the lifetime is given by 

r = h/(Tc 2 ) (49) 
If the particle a has a definite non-zero momentum p, then 

V'(q) = ^(q-p) (so) 

and 

oo 

W| p) (0,t) = I J dm\^{m)\ 2 e'^ t \ 2 (51) 

In a number of works [H EJ IE] it was noticed that if one interprets the state |p) 
(where |p| = macsinh^) as a state of unstable particle moving with definite speed 
|v| = ctanh#, then the decay law of the moving particle ()51)1 cannot be connected 
with the non-decay law of the particle at rest by Einstein's formula Q, i.e., 

W| P) (0,*) ^ cU| O) (O,t/cosh0) (52) 

This observation prompted authors of jHEHEj to question the applicability of Einstein's 
special relativity to particle decays. However, at a closer inspection it appears that this 
result does not contradict Einstein's time dilation formula directly. Formula 
refers to observations made on the same particle from two frames of reference moving 
with respect to each other. If from the point of view of observer O the particle is 
described by the state vector |0) which has zero momentum and zero velocity, then 
from the point of view of O' this particle is described by the state 

|p) = e ^ Ke |0) (53) 

16 By noting that for the particle at rest its energy is identified with mc 2 , this formula can be 
compared, for example, with eq. (3.8) in 
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which is an eigenstate of the velocity operator 



V = c 2 P /H (54) 

but is not an eigenstate of the momentum operator P . Therefore, strictly speaking, 
its non-decay law is not described by u\ p )(0,t). In order to compare with Einstein's 
formula (jUJ), we need to calculate the decay law in the moving frame of reference 
uj\y}(9,t). This is done in section HU 



10 Decays caused by boosts 

Let us now discuss the non-decay probability u\^,)(6, 0) at t = in the moving frame of 
reference. Instead of eq. (147(1 . it is more convenient to use the general definition (|22J1 
which expresses {6, 0) as the square of the norm of the projection of \^f(9, 0)) on the 
subspace 7i a . We are going to prove that for a normalized this probability cannot 
be equal to 1 for all non-zero 9. Suppose that this statement is wrong, so that (in 
Agreement with Assumption V) for any £ 7i a and any 9 > 0, the vector e^ Kz9 \^) 
belongs to 7i a . Then the subspace H a is invariant under action of boosts e^ Kz0 and 
operator K z commutes with the projection T. Then from the Poincare commutator 
(fT5|) and [T, (Pq) z ] = it follows by Jacobi identity that 



= j-[K z ,[T,(P ) z }]- j-[(P ) z ,[T,K z }} 
= 

which contradicts eq. (|19j). This contradiction implies that the state e^ Kz8 \^l) does 
not correspond to the particle a with 100% probability; this state must contain con- 
tributions from the decay products even at t = 

17 Indeed, taking into account V z \0) « and eqs. (jT7|l - i(T%|) . we obtain 



V z e* K >° 



0) 



e e " V z e h 



|0) 



= P H 



4-A" = 



,v t 



- c tanh 9 

V^tanhff I") 



: tanh ( 



2 K z t 



|0) 
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e T?*-»|tf) £ H a (55) 
w w (0,O) < 1, for 9^0 (56) 

This is the "decay caused by boost" . Thus we found that Assumption V from section HI 
is not correct. Sometimes this assumption is formulated as [13] "Flavor is the quantum 
number that distinguishes the different types of quarks and leptons. It is a Lorentz 
invariant quantity. For example, an electron is seen as an electron by any observer, 
never as a muon. " Although the above statement about the electron is correct (because 
the electron is a stable particle), this is not true about the muon which, according to 
(I55j) can be seen as a single particle by the observer at rest and as a group of three 
decay products (an electron, a neutrino, and an antineutrino) by the moving observer. 



11 Decay law in the moving reference frame 

Unfortunately exact evaluation of eq. (jTFj) for 9 ^ is not possible, so we need to make 
approximations. Let us discuss properties of the initial state |\&) in more detail. First, 
in all realistic cases this state is not an eigenstate of the total momentum operator, so 
the wave function is not localized at one point in the momentum space (as assumed 
in fl5()j) ) but has a spread (or uncertainty) of momentum |Ap| and, correspondingly an 
uncertainty of position |Ar| ~ ft/|Ap|. On the other hand, the state G 7i a is not 
an eigenstate of the mass operator M, and is characterized by the uncertainty of 
mass T (see Fig. 1) that is related to the lifetime of the particle r by formula ft 
is important to note that in all cases of practical interest the mentioned quantities are 
related by inequalities 



|Ap| > Tc (57) 
|Ar| < cr (58) 

In particular, the latter inequality means that the uncertainty of position is mush less 
than the distance passed by light during the lifetime of the particle. For example, 
in the case of muon r ~ 2.2 • lCT 6 s and, according to (fo'Hj). the spread of the wave 
function in the position space must be much less than 600m, which is a reasonable 
assumption. Therefore, we can safely assume that the factor |/x(m)| 2 in (|^7j) has a 
sharp peak near the value m = m a . Then we can move the value of the smooth 
function y^^/>(Aq)7(Aq, m)7*(q, m) at m = m a outside the integral by m. 18 

18 Note that if we assumed instead of l(57l) that the spread of the momentum-space wave function is 
much less than Tc we could have moved the value of the smooth function |/x(mo)| 2 outside the integral 



22 



u m (9,t) « | C /qU^(Aq) 7 (Aq,m a ) 7 *(q,m a )| 2 | J dm|/i(m)| 2 e -^f 

q m 6 +m c 

= J rfq-^|^(Lq)| 2 | J rfm|/i(m)|V^f 

q m b +m c 

CO 

= J dp\^{p)\ 2 \ J dmlfiim^e-^L-^l 2 (59) 

m b +m c 

Now we are going to use formula (jHHj) to calculate the decay law for an approximate 
eigenstate of momentum in 7i a whose wave function ip(p) is localized near zero mo- 
mentum p = 0, though inequality (|57|) still remains valid. We denote this state vector 
by symbol |0]. The state ef K * e \0] is an approximate eigenstate of the velocity opera- 
tor (J53|) for all values of 6. 19 Therefore U\ ](9,t) in eq. (|50|l can be understood as the 
non-decay law of a particle with definite speed v = ctanh^. 

Note that the second factor in the integrand in is a slowly varying function of 
p. Therefore, we can set |^(p)| 2 ~ S(p) in eq. (|59|) and obtain 



00 

u\ 0] (9,t) w I J dm\fi{m)\ 2 e-^^ m2c4+m ' c4sinh2e \ 2 (60) 

mh+m c 

If we approximately 20 identify m a c sinh# with the momentum |p| of the particle a from 
the point of view of the moving observer O' then 



00 

u m (9,t) « I J dm\fi(m)\ 2 e-^*\ 3 (61) 

So, in this approximation the non-decay law ()61|) in the frame of reference O' moving 
with the speed ctanh^ takes the same form as the non-decay law ()51|) of a particle 
moving with momentum m a c sinh# with respect to the stationary observer O. 

by m. This would mean that the non-decay law in the moving frame of reference is controlled by 
the spread of the particle wavefunction in the momentum space rather than by its mass uncertainty 
which disagrees with experimental observations. 

19 See footnote 17. 

20 See remark after eq. (|5*5|l . 
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12 Numerical results 



In this section we are going to perform numerical calculations of the differences between 
the classical Einstein's formula (JHJ) 

^HM = -|o](0,^) (62) 

and the actual non-decay law (ffiTij) of a moving particle having either definite mo- 
mentum or sharply defined speed. In these calculations we assumed that the mass 
distribution |/i(m)| 2 of the unstable system (see eqs. (J4*3j) - (}4"4"j) and Fig. 1) is centered 
at the value of mass m a = 1000 MeV/c 2 , the total mass of the decay products was 
nib + f^c = 900 MeV/c 2 , and the width of the mass distribution was T= 20 MeV/c 2 . 
These values do not correspond to any real particle, but they are typical for strongly 
decaying baryons. 

It is convenient to measure time in units of the classical lifetime Tq cosh^. Denoting 
X = t/(T cosh#), we find that Einsteinian non-decay laws (p32*|) for any value of 9 are 
given by the same universal function u dass (x)- This function was calculated for values 
of x m the interval from to 6 with the step of 0.1 The calculations were performed by 
direct numerical integration of eq. (f4"KJ) using the Mathematica program shown below 

gamma = 20 
mass = 1000 
theta =0.0 

Do [Print [(1/0. 9375349) Abs [NIntegrate [gamma/ (2 Pi) / (gamma~2/4 

+(x - mass) "2) Exp[ I t Sqrt [x~2 + mass~2 (Sinn [theta]) "2] 

Cosh [theta] / gamma], {x, 900, 1010, 1100, 300000}, MinRecursion -> 3, 

MaxRecursion -> 16, PrecisionGoal -> 8, WorkingPrecision -> 18]] "2], 

{t, 0.0, 6.0, 0.1}] 

As expected, function u class (x) is ver Y close to the exponent e~ x . 21 Next we used eq. 
(JbTJ)) and the above Mathematica program to calculate the non-decay laws a>| ] (0, x) 
in moving reference frames for three values of the parameter 9 (=theta), namely 
0.2, 1.4, and 10.0, that correspond to velocities of 0.197c, 0.885c, and 0.999999995c, 
respectively. These calculations qualitatively confirmed the time dilation formula ()(i2)l 

21 This function is represented by a thick solid line in Fig. 2. The magnitudes of uj class (x) for small 
values of the argument % are too large to be shown on the scale of Fig. 2. 
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Figure 2: Corrections to the Einstein's "time dilation" formula for the non-decay 
law of unstable particle moving with the speed v = ctanh#. Parameter x is time 
measured in units of To/ cosh 9. 

to the accuracy of better than 0.3%. However, they also revealed important differences 
w\o]{9, x) ~ uj class (x) which are plotted as thin lines in Fig. 2. 22 

The lifetime of the particle a considered in our example (r ~ 2 x 10 -22 s) is too 
short to be observed experimentally. So, calculated corrections to the Einstein's time 
dilation law have only illustrative value. However, from these data we can estimate 
the magnitude of corrections for particles whose time-dependent non-decay laws can 
be measured in a laboratory, e.g., for muons (r pa 2 x lCT 6 s, r ~ 2 x 10~ 9 eV/c 2 , 
m a pa 105MeV/c 2 ). Taking into account that the magnitude of corrections is roughly 
proportional to the ratio T/m a [H Ej we can expect that for muons the maximum 
correction should be about 2 x 10 -18 which is much smaller than the accuracy of 
modern experiments (about 10~ 3 |2j). So, experimental observation of the predicted 
corrections requires significant improvements of existing experimental techniques. 

Note that formulas ()60|) and (JoTj) for the particle seen from the moving frame are 
approximate in the sense that they ignore the "decay caused by boost" expressed by 
eq. (jHEJ). Nevertheless, our major approximation (|H7|) is well justified and it cannot 
explain the discrepancy of our results (Einsteinian time dilation of the decay plus small 
interaction-dependent corrections) with the conclusion of ref. jTj about the acceleration 
of the decay of moving particles. In our view this conclusion does not refer to the 
experimentally measured non-decay law that is defined as the probability of finding 

22 Results presented in Fig. 2 are different from those in Fig. 1 in ref. ^ due to the more accurate 
calculation procedure employed in the present work. 
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one unstable particle. Instead of eq. ([21)1. the non-decay probability was defined in [7j 
by formula 



ui\m){9, t) 



K*(M)|e 



-z m \V(9, 0))| 2 



whose physical meaning remains unclear. 

13 Particle decays in different Dirac's forms of dy- 
namics 

In this article we assumed that the interaction responsible for the particle decay has the 
Bakamjian-Thomas instant form of dynamics. This assumption was used to simplify 
calculations, but there is no good reason to believe that real interactions have this 
form. Then, naturally, one may ask a question: "Is there another form of dynamics 
in which Einstein's time dilation formula is exactly true?" The answer is No. 
In any instant form of dynamics (including no n-Bakamjian- Thomas instant forms of 
dynamics) the boost operators contain interaction terms, so the results (jo3j) - (p)6r) are 
still valid, and eq. (jHI) holds only approximately. In the point form dynamics 23 the 
subspace H a of the unstable particle is invariant with respect to boosts, so there can 
be no boost-induced decays (fo3j). However, due to the interaction-dependence of the 
total momentum operator P, one should expect decays induced by space translations 



This prediction is not confirmed by experiments which show that the composition of 
an unstable particle is not affected by space translations of the observer. 24 Therefore, 
the point form of dynamics is not acceptable for the description of decays. Similarly, 
never observed translation- and/or rotation-induced decays are characteristic for all 
non- instant forms of dynamics, e.g., the front form. Therefore only the instant form of 
dynamics is appropriate for the description of particle decays. Note, however, that the 
decay slowdown and the absence of translation-induced decays are firmly established 
only for unstable systems whose lifetimes are long enough to be measured experi- 
mentally, i.e., those governed by electromagnetic and weak interactions. For strongly 
interacting resonances these properties are beyond experimental resolution, and the 
non- instant forms of dynamics cannot be ruled out. 

23 where generators of boosts Ko and rotations Jo are kinematical, while generators of space trans- 
lations P and time translations H contain interaction terms. 

24 An even more striking contradiction between predictions of the point form dynamics and observa- 
tions is that the decay of moving particles accelerates cosh 9 times instead of experimentally observed 
slowdown 




H a , for a^O 



(63) 
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14 Summary 



In this paper we analyzed the relationships between the non-decay laws in the moving 
reference frame and in the reference frame at rest. We used a rigorous quantum 
relativistic approach that is applicable to any unstable system independent on the 
nature of interaction governing the decay. A complete description of dynamics in 
different reference frames was obtained by using relativistic Postulates I - VI and rules 
of quantum mechanics only. We found that Assumptions V and VI of special relativity 
are not needed. Moreover, their consequences (the universal and exact time dilation of 
the decay of moving particles) are in contradiction with rigorous calculations. Although 
the time dilation (jHJ) of special relativity was qualitatively confirmed by our results, we 
also found small corrections to this formula that depend on the strength of interaction. 
In a broader sense our results indicate that clocks viewed from the moving reference 
frame do not go exactly cosh 9 slower. The exact amount of time dilation depends on 
the physical makeup of the clock and on interactions responsible for the operation of 
the clock. 

I am grateful to Dr. M. I. Shirokov for reading the manuscript, valuable comments, 
and discussions. 
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